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Abstract 

Grand Unified Theories predict relationships between the GUT-scale quark and lepton masses. Using 
new data in the context of the MSSM, we update the values and uncertainties of the masses and mixing 
angles for the three generations at the GUT scale. We also update fits to hierarchical patterns in the 
GUT-scale Yukawa matrices. The new data shows not all the classic GUT-scale mass relationships remain 
in quantitative agreement at small to moderate tan /3. However, at large tan /3, these discrepancies can be 
eliminated by finite, tan /3-enhanced, radiative, threshold corrections if the gluino mass has the opposite 
sign to the wino mass. 

Explaining the origin of fermion masses and mixings remains one of the most important goals in our 
attempts to go beyond the Standard Model. In this, one very promising possibility is that there is an 
underlying stage of unification relating the couplings responsible for the fermion masses. However we are 
hindered by the fact that the measured masses and mixings do not directly give the structure of the underlying 
Lagrangian both because the data is insufficient unambiguously to reconstruct the full fermion mass matrices 
and because radiative corrections can obscure the underlying structure. In this letter we will address both 
these points in the context of the MSSM. 

We first present an analysis of the measured mass and mixing angles continued to the GUT scale. The 
analysis updates previous work, using the precise measurements of fermion masses and mixing angles from 
the b-factories and the updated top-quark mass from CDF and DO. The resulting data at the GUT scale 
allows us to look for underlying patterns which may suggest a unified origin. We also explore the sensitivity 
of these patterns to tan /3-enhanced, radiative threshold corrections. 

We next proceed to extract the underlying Yukawa coupling matrices for the quarks and leptons. There 
are two difficulties in this. The first is that the data cannot, without some assumptions, determine all elements 
of these matrices. The second is that the Yukawa coupling matrices are basis dependent. We choose to work 
in a basis in which the mass matrices are hierarchical in structure with the off-diagonal elements small relative 
to the appropriate combinations of on-diagonal matrix elements. This is the basis we think is most likely to 
display the structure of the underlying theory, for example that of a spontaneously broken family symmetry 
in which the hierarchical structure is ordered by the (small) order parameter breaking the symmetry. With 
this structure to leading order the observed masses and mixing angles determine the mass matrix elements 
on and above the diagonal, and our analysis determines these entries, again allowing for significant tan/3 
enhanced radiative corrections. The resulting form of the mass matrices provides the "data" for developing 
models of fermion masses such as those based on a broken family symmetry. 

The data set used is summarized in Table [TJ Since the fit of reference [4 (RRRV) to the Yukawa texture 
was done, the measurement of the Standard-Model parameters has improved considerably. We highlight 
a few of the changes in the data since 2000: The top-quark mass has gone from M t = 174.3 ± 5 GeV to 
M t = 170.9 ± 1.9 GeV. In 2000 the Particle Data Book reported m b (m b ) = 4.2 ± 0.2 GeV [5] which has 
improved to mb(mi,) = 4.2 ± 0.07 GeV today. In addition each higher order QCD correction pushes down 
the value of m^iMz) at the scale of the Z bosons mass. In 1998 m b (Mz) = 3.0 ± 0.2 GeV [6] and today 
it is mt,(Mz) = 2.87 ± 0.06 GeV [7]. The most significant shift in the data relevant to the RRRV fit is a 
downward revision to the strange-quark mass at the scale hl = 2 GeV from iti s ({1l) w 120 ± 50 MeV [5] to 
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Low-Energy Parameter 


Value(Uncertainty in last digit (s)) 


Notes and Reference 




U. 40(10 ) 


y ij\ ) ri/Siniiaiion ij 


m s (fi L )/md(fJ.L) 


19.0(1.0 ) 


FUri HiStimation [ij 


m u {iiL) + m d (fi L ) 


|o.o(o.uj, f.o(i.ojj iviev 


r^Uri, l^uaiK Masses, pg lo 






[IJ. [ INon-lattice, Lattice ) 


V m d~~ m u 


22.8(4) 


Martemyanov and Sopov [5] 


m s (fi L ) 


[103(20), 95(20)] MeV 


PDB, Quark Masses, pg 15 






[1]. [Non-lattice, lattice] 


m u (fi L ) 


3(1) MeV 


PDB, Quark Masses, pg 15 






PP. Non-lattice. 


m d (i± L ) 


6.0(1.5) MeV 


PDB, Quark Masses, pg 15 






pp. Non-lattice. 


m c (m c ) 


1.24(09) GeV 


PDB, Quark Masses, pg 16 






PP. Non-lattice. 


m b (m b ) 


4.20(07) GeV 


PDB, Quark Masses, pg 16,19 






PP. Non-lattice. 


M f 


170.9 (1.9) GeV 


CDF & DO Pole Mass 


(M e ,M M M T ) 


(0.511(15), 105.6(3.1), 1777(53) ) MeV 


3% uncertainty from neglect- 






ing Y e thresholds. 


A Wolfenstein parameter 


0.818(17) 


PDB Ch 11 Eq. 11.25 pp 


p Wolfenstein parameter 


0.221(64) 


PDB Ch 11 Eq. 11.25 [T] 


A Wolfenstein parameter 


0.2272(10) 


PDB Ch 11 Eq. 11.25 pp 


rj Wolfenstein parameter 


0.340(45) 


PDB Ch 11 Eq. 11.25 [T] 




/0.97383(24) 0.2272(10) 0.00396(09) \ 




\Vckm\ 


0.2271(10) 0.97296(24) 0.04221(80) 


PDB Ch 11 Eq. 11.26 pp 




\0.00814(64) 0.04161(78) 0.999100(34)/ 




sin 2(3 from CKM 


0.687(32) 


PDB Ch 11 Eq. 11.19 pp 


Jarlskog Invariant 


3.08(18) x 10~ 5 


PDB Ch 11 Eq. 11.26 pp 


VHiggs(Mz) 

( d£ M {M z ), a s (M z ), sm 2 e w {M z ) ) 


246.221(20) GeV 
( 127.904(19), 0.1216(17), 0.23122(15)) 


Uncertainty expanded, pp 
PDB Sec 10.6 pp 



Table 1: Low-energy observables. Masses in lower-case m are MS running masses. Capital M indicates pole 
mass. The light quark's (u,d,s) mass are specified at a scale hl = 2 GeV. Vckm are the Standard Model's 
best fit values. 



today's value m s (/ii) = 103 ± 20 MeV. We also know the CKM unitarity triangle parameters better today 
than six years ago. For example, in 2000 the Particle Data book reported sin 2(3 = 0.79 ± 0.4 [5] which is 
improved to sin 2(3 — 0.69 ± 0.032 in 2006 1 . The sin 2/3 value is about 1.2 a off from a global fit to all the 
CKM data [H], our fits generally lock onto the global-fit data and exhibit a 1 a tension for sin 2(3. Together, 
the improved CKM matrix observations add stronger constraints to the textures compared to data from 
several years ago. 

We first consider the determination of the fundamental mass parameters at the GUT scale in order simply 
to compare to GUT predictions. The starting point for the light-quark masses at low scale is given by the 
X 2 fit to the data of Table CD 

m u (fJ, L ) = 2.7 ±0.5 MeV m d (fi L ) = 5.3 ± 0.5 MeV m s (fi L ) = 103 ± 12 MeV. (1) 

Using these as input we determine the values of the mass parameters at the GUT scale for various choices of 
tan/3 but not including possible tan/3 enhanced threshold corrections. We do this using numerical solutions 
to the RG equations. The one-loop and two-loop RG equations for the gauge couplings and the Yukawa 
couplings in the Standard Model and in the MSSM that we use in this study come from a number of sources 
[5] [S] [TU] [H] • The results are given in the first five columns of Table [5J These can readily be compared to 
expectations in various Grand Unified models. The classic prediction of SU(5) with third generation down- 
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quark and charged-lepton masses given by the coupling B 5f.l0f.5n El is rnb(Mx)/m T (Mx) — 1 [II]- This 
ratio is given in Table [2] where it may be seen that the value agrees at a special low tan/3 value but for large 
tan j3 it is some 25% smaller than the GUT prediction^]. A similar relation between the strange quark and 
the muon is untenable and to describe the masses consistently in SU (5) Georgi and Jarlskog [ 14] proposed 
that the second generation masses should come instead from the coupling C 5/. 10/. A5h leading instead to 
the relation 3m s (Mx) /m^Mx) — 1. As may be seen from Table 2 in all cases this ratio is approximately 
0.69(8). The prediction of Georgi and Jarlskog for the lightest generation masses follows from the relation 
Det{M d ) / Det{M l ) = 1. This results from the form of their mass matrix which is given bjH 

A' 

M d = I A C I , M l = I A -3C I (2) 




in which there is a (1, 1) texture zero3 and the determinant is given by the product of the (3, 3), (1, 2) and 
(2, 1) elements. If the (1, 2) and (2, 1) elements are also given by 5/.10/.5# couplings they will be the same 
in the down-quark and charged-lepton mass matrices giving rise to the equality of the determinants. The 
form of eq((2]) may be arranged by imposing additional continuous or discrete symmetries. One may see from 
Table [2] that the actual value of the ratio of the determinants is quite far from unity disagreeing with the 
Georgi Jarlskog relation. 

In summary the latest data on fermion masses, while qualitatively in agreement with the simple GUT 
relations, has significant quantitative discrepancies. However the analysis has not, so far, included the SUSY 
threshold corrections which substantially affect the GUT mass relations at large tan /3 [15] . A catalog of the 
full SUSY threshold corrections is given in [16] . The particular finite SUSY thresholds discussed in this letter 
do not decouple as the super partners become massive. We follow the approximation described in Blazek, 
Raby, and Pokorski (BRP) for threshold corrections to the CKM elements and down-like mass eigenstates 
|17j . The finite threshold corrections to Y e and Y u and are generally about 3% or smaller 

SY U , SY d < 0.03 (3) 

and will be neglected in our study. The logarithmic threshold corrections are approximated by using the 
Standard-Model RG equations from Mz to an effective SUSY scale Mg. 

The finite, tan /3-enhanced Y d SUSY threshold corrections are dominated by the a sbottom-gluino loop, 
a stop-higgsino loop, and a stop-chargino loop. Integrating out the SUSY particles at a scale Ms leaves the 
matching condition at that scale for the Standard-Model Yukawa couplings: 

Sm sch Y uSM = sin/3 Y u (4) 
Sm sch Y d SM = cos E/f ( 1 + r d + 4 M T u Vckm) Y d U d R (5) 



CKM *CKM J -'diag^fl 

Y eSM = cos/3 Y e . (6) 

All the parameters on the right-hand side take on their MSSM values in the DR scheme. The factor 8m sc h 
converts the quark running masses from MS to DR scheme. The j3 corresponds to the ratio of the two 
Higgs VEVs v u /vd — tan/3. The U matrices decompose the MSSM Yukawa couplings at the scale Ms: 
yu = U^Y£ ag U% and Y d = U^Y^U^. The matrices Y£ ag and Y& ag are diagonal and correspond to 
the mass eigenstates divided by the appropriate VEV at the scale Ms- The CKM matrix is given by 
Vckm — U L U L \ The left-hand side involves the Standard-Model Yukawa couplings. The matrices T w and 
T d encode the SUSY threshold corrections. 

If the squarks are diagonalized in flavor space by the same rotations that diagonalize the quarks, the 
matrices T u and T d are diagonal: T d = diag(7d, 7<j, 7h), T u = diag(7„, j u , 7 t ). In general the squarks are 



1 5j-, 10 f refer to the SU(5) representations making up a family of quarks and leptons while 5h is a five dimensional 
representation of Higgs scalars. 

2 We'd like to thank Ilja Dorsner for pointing out that the tan (3 dependence of mi 1 /m T (Mx) is more fiat than in previous 
studies (e.g. ref. [13]). This change is mostly due to the higher effective SUSY scale Mg, the higher value of ce s (Mz) found in 
global standard model fits, and smaller top-quark mass Mt- 

3 The remaining mass matrix elements may be non-zero provided they do not contribute significantly to the deteminant 

4 Below we discuss an independent reason for having a (1, 1) texture zero. 
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Parameters 


Input SUSY Parameters 


tan j3 


1.3 


10 


38 


50 


38 


38 


lb 














-0.22 


+0.22 


Id 














-0.21 


+0.21 


It 

















-0.44 


Parameters 


Corresponding GUT-Scale Parameters with Propagated Uncertainty 


y^Mx) 
y h (M x ) 

y T (M x ) 

(m u /m c ){M x ) 
{m d /m s ){Mx) 
{m e /m^){M x ) 
(m c /m t )(M x ) 
(m s /m b )(M x ) 
(m M /m T )(M x ) 

A(M X ) 

\{M X ) 

p{M x ) 

fj(Mx) 
J{M X ) x 10 


0.0113™ 2 
0.0114(3) 
0.0027(6) 
0.051(7) 
0.0048(2) 

00 09+ 001 

u ' uuuy -0. 00006 

0.014(4) 
0.059(2) 

U - OD -0.01 

0.227(1) 
0.22(6) 
0.33(4) 

1 4+2-2 


0.48(2) 
0.051(2) 
0.070(3) 
0.0027(6) 
0.051(7) 
0.0048(2) 
0.0025(2) 
0.019(2) 
0.059(2) 
0.77(2) 
0.227(1) 
0.22(6) 
0.33(4) 
2.6(4) 


0.49(2) 
0.23(1) 
0.32(2) 
0.0027(6) 
0.051(7) 
0.0048(2) 
0.0024(2) 
0.017(2) 
0.054(2) 
0.75(2) 
0.227(1) 
0.22(6) 
0.33(4) 
2.5(4) 


0.51(3) 
0.37(2) 
0.51(4) 
0.0027(6) 
0.051(7) 
0.0048(2) 
0.0023(2) 
0.016(2) 
0.050(2) 
0.72(2) 
0.227(1) 
0.22(6) 
0.33(4) 
2.3(4) 


0.51(2) 
0.34(3) 
0.34(2) 
0.0026(6) 
0.051(7) 
0.0048(2) 
0.0023(2) 
0.018(2) 
0.054(2) 
0.73(3) 
0.227(1) 
0.22(6) 
0.33(4) 
2.3(4) 


0.51(2) 
0.34(3) 
0.34(2) 
0.0026(6) 
0.051(7) 
0.0048(2) 
0.0023(2) 
0.010(2) 
0.054(2) 
0.46(3) 
0.227(1) 
0.22(6) 
0.33(4) 
1.0(2) 


Parameters 


Comparison with GUT Mass Ratios 


(m b /m T )(M x ) 
(3m s /m M )(M x ) 
(m d /3m e )(M x ) 


i.oolH 4 

0.82(7) 


0.73(3) 
0.69(8) 
0.83(7) 


0.73(3) 
0.69(8) 
0.83(7) 


0.73(4) 
0.69(8) 
0.83(7) 


1.00(4) 
0.9(1) 
1.05(8) 


1.00(4) 
0.6(1) 
0.68(6) 




^7+ 08 

u -°'-0.26 


0.42(7) 


0.42(7) 


0.42(7) 


0.92(14) 


0.39(7) 



Table 2: The mass parameters continued to the GUT-scale Mx for various values of tan/3 and threshold 
corrections 74,6^. These are calculated with the 2-loop gauge coupling and 2-loop Yukawa coupling RG 
equations assuming an effective SUSY scale Ms = 500 GeV. 



not diagonalized by the same rotations as the quarks but provided the relative mixing angles are reasonably 
small the corrections to flavour conserving masses, which are our primary concern here, will be second order 
in these mixing angles. We will assume T u and T d are diagonal in what follows. 

Approximations for r u and r d based on the mass insertion approximation are found in [TB] [Tj5] [5U] : 

2 a t tan/? . 9 2 2\ 2 tan /3 \iA^ . . 

7 * ~ Vt M I^2- /3 K+ m tV^ ) ~ Vt oo_2 ^2 ( 7 ) 



7« ~ -9lM 2 p^^l3(ml 1 ,m 2 X2 ,ml) ~ (8) 



8 2 tan/3 2 2 „ # 2% 4 2 ta - n Pv M 3 

3 167r2 °i °2 3 167r2 m- 

b 

8 2 tan£ 2 2 , f 2N 4 2 tan / ? M^3 

3 1671-2 fl i "2 3 167r2 m~ 



where I3 is given by 



, o o 9n a 2 fe 2 log w + & 2 c 2 log t- + c 2 a 2 log % 
^'•^ - ( a 2 - b 2 )(6 2 - c 2 K a 2 - C 2 ) " ■ (U) 



In these expressions q refers to superpartner of q. x* indicate chargino mass eigenstates. n is the coefficient 
to the H u H d interaction in the superpotential. M 1; M 2 , M 3 are the gaugino soft breaking terms. A* refers to 
the soft top-quark trilinear coupling. The mass insertion approximation breaks down if there is large mixing 
between the mass eigenstates of the stop or the sbottom. The right-most expressions in eqs (l7l9U0[) assume 
the relevant squark mass eigenstates are nearly degenerate and heavier than M3 and /1. These expressions ( 
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eqs [7] - 1101) provide an approximate mapping from a supersymmetric spectra to the 7$ parameters through 
which we parameterize the threshold corrections; however, with the exception of Column A of Table 2] we 
do not specify a SUSY spectra but directly parameterize the thresholds corrections through 7,. 

The separation between 7^ and jd is set by the lack of degeneracy of the down-like squarks. If the 
squark masses for the first two generations are not degenerate, then there will be a corresponding separation 
between the (1,1) and (2,2) entries of T d and T u . If the sparticle spectra is designed to have a large A* and 
a light stop, 7t can be enhanced and dominate over jf,. Because the charm Yukawa coupling is so small, the 
scharm-higgsino loop is negligible, and "f u follows from a chargino squark loop and is also generally small 
with values around 0.02 because of the smaller g2 coupling. In our work, we approximate TVfa ~ ~ 0. 
The only substantial correction to the first and second generations is given by 7^ |15j . 

As described in BRP, the threshold corrections leave \V US \ and IV^/V^I unchanged to a good approxi- 
mation. Threshold corrections in T u do affect the V u b and V c b at the scale M$ giving 

ySM _ yMSSM ySM _ yMSSM 

v ub v ub ^ cb • cb ^ / \ (~\0\ 

yMSSM — yMSSM ~ W* ' u > ' 

v ub v cb 

The threshold corrections for the down-quark masses are given approximately by 



m d 
m. 



m°(l + 7d + 7n)~ 1 



nib — m b (1 + jb + It) 

where the superscript denotes the mass without threshold corrections. Not shown are the nonlinear effects 
which arise through the RG equations when the bottom Yukawa coupling is changed by threshold effects. 
These are properly included in our final results obtained by numerically solving the RG equations. 

Due to our assumption that the squark masses for the first two generations are degenerate, the combina- 
tion of the GUT relations given by (det M l / det M d ) (3 m s /m p ) 2 (rrib/m T ) = 1 is unaffected up to nonlinear 
effects. Thus we cannot simultaneously fit all three GUT relations through the threshold corrections. A best 
fit requires the threshold effects given by 

7b + 7f ps -0.22 ±0.02 (13) 
ld + lu « -0.21 ±0.02. (14) 

giving the results shown in the penultimate column of Table[2j just consistent with the GUT predictions. The 
question is whether these threshold effects are of a reasonable magnitude and, if so, what are the implications 
for the SUSY spectra which determine the 7^? From eqs (|9|10jl . at tan/3 = 38 we have 



2 

b 



— — - -0.5, — §• - 1.0 

mf m~ 

b d 

The current observation of the muon's (g — 2) M is 3.4 a [21] away from the Standard-Model prediction. 
If SUSY is to explain the observed deviation, one needs tan/3 > 8 [22] and \iMi > [23] , With this sign 
we must have negative and the d, Is squarks only lightly split from the b squarks. M3 negative is 

characteristic of anomaly mediated SUSY breaking[24 and is discussed in [25 t 26j [20J [27 . Although we have 
deduced M3 < from the approximate eas (|9U0p . the the correlation persists in the near exact expression 
found in eq(23) of ref [17) . Adjusting to different squark splitting can occur in various schemes 28 . However 
the squark splitting can readily be adjusted without spoiling the fit because, up to nonlinear effects, the 
solution only requires the constraints implied by eq (|13p . so we may make 7^ > jd and hence make m? < m| 
by allowing for a small positive value for 7 t . In this case A* must be positive. 

It is of interest also to consider the threshold effects in the case that //M3 is positive. This is illustrated 
in the last column of Table [2] in which we have reversed the sign of jd, consistent with positive [iM^ , and 
chosen 7^ ~ 7^ as is expected for similar down squark masses. The value of 74 is chosen to keep the equality 
between nib and m T . One may see that the other GUT relations are not satisfied, being driven further away 
by the threshold corrections. Reducing the magnitude of 7?, and 7^ reduces the discrepancy somewhat but 
still limited by the deviation found in the no-threshold case (the fourth column of Table \2§ ■ 
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Parameter 


2001 RRRV 


Fit AO 


Fit BO 


Fit Al 


Fit Bl 


Fit A2 


Fit B2 


tan/3 


Small 


1.3 


1.3 


38 


38 


38 


38 


a! 


0(1) 














-2.0 


-2.0 


e u 


0.05 


0.030(1) 


0.030(1) 


0.0491(16) 


0.0491(15) 


0.0493(16) 


0.0493(14) 


td 


0.15(1) 


0.117(4) 


0.117(4) 


0.134(7) 


0.134(7) 


0.132(7) 


0.132(7) 


\b'\ 


1.0 


1.75(20) 


1.75(21) 


1.05(12) 


1.05(13) 


1.04(12) 


1.04(13) 


arg(b') 


90° 


+ 93(16)° 


-93(13)° 


+ 91(16)° 


-91(13)° 


+ 93(16)° 


-93(13)° 


a 


1.31(14) 


2.05(14) 


2.05(14) 


2.16(23) 


2.16(24) 


1.92(21) 


1.92(22) 


b 


1.50(10) 


1.92(14) 


1.92(15) 


1.66(13) 


1.66(13) 


1.70(13) 


1.70(13) 


\c\ 


0.40(2) 


0.85(13) 


2.30(20) 


0.78(15) 


2.12(36) 


0.83(17) 


2.19(38) 


arg(c) 


-24(3)° 


-39(18)° 


-61(14)° 


-43(14)° 


-59(13)° 


-37(25)° 


-60(13)° 



Table 3: Results of a \ 2 fit of eqs (|15|16|) to to the data in Table [2] in the absence of threshold corrections. 
We set a' as indicated and set d = d! = d = and / = /' = 1 at fixed values. 



At tan/3 near 50 the non-linear effects are large and b — t unification requires 7b + 7t ~ —0.1 to —0.15. 
In this case it is possible to have t — b — r unification of the Yukawa couplings. For fi > 0, M3 > 0, the 
"Just-so" Split-Higgs solution of references [29J, [30j [HI [32] can achieve this while satisfying both b — > s 7 and 
(g — 2)^ constraints but only with large 7b and jt and a large cancellation in 7b + 7*. In this case, as in the 
example given above, the threshold corrections drive the masses further from the mass relations for the first 
and second generations because fj, M3 > 0. It is possble to have t — b — r unification with /x M3 < 0, satisfying 
the b — > s 7 and (g — 2) M constraints in which the GUT predictions for the first and second generation of 
quarks is acceptable. Examples include Non-Universal Gaugino Mediation [33] and AMSB; both have some 
very heavy sparticle masses ( > 4 TeV) [20] . Minimal AMSB with a light sparticle spectra( < 1 TeV) , while 
satisfying [g — 2)^ and b — ► s 7 constraints, requires tan/3 less than about 30 [23] . 

We turn now to the second part of our study in which we update previous fits to the Yukawa matrices 
responsible for quark and lepton masses. As discussed above we choose to work in a basis in which the 
mass matrices are hierarchical with the off-diagonal elements small relative to the appropriate combinations 
of on-diagonal matrix elements. This is the basis we think is most likely to display the structure of the 
underlying theory, for example that of a spontaneously broken family symmetry, in which the hierarchical 
structure is ordered by the (small) order parameter breaking the symmetry. With this structure to leading 
order in the ratio of light to heavy quarks the observed masses and mixing angles determine the mass matrix 
elements on and above the diagonal provided the elements below the diagonal are not anomalously large. 
This is the case for matrices that are nearly symmetrical or for nearly Hermitian as is the case in models 
based on an 5O(10) GUT. 

For convenience we fit to symmetric Yukawa coupling matrices but, as stressed above, this is not a critical 
assumption as the data is insensitive to the off-diagonal elements below the diagonal and the quality of the 
fit is not changed if, for example, we use Hermitian forms. We parameterize a set of general, symmetric 
Yukawa matrices as: 

Y U (M X ) = y^lb'el f 4 a'e| I , (15) 



d'et 


b'el 


c'el 


b'el 


f'el 


a' el 


c'el 


a' el 


1 


d4 


be d 


ce d \ 




fe\ 


.4 


c4 


a e\ 





Y d (M x ) = 2/33 I be d fej aej\. (16) 
Ve d ael l] 

Although not shown, we always choose lepton Yukawa couplings at Mx consistent with the low-energy lepton 
masses. Notice that the / coefficient and e d are redundant (likewise in Y u ). We include / to be able to 
discuss the phase of the (2,2) term. We write all the entries in terms of e so that our coefficients will be 
0(1). We will always select our best e parameters such that |/| = 1. 

RRRV noted that all solutions, to leading order in the small expansion parameters, only depend on two 
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Parameter 


A 


B 


C 


B2 


C2 




tan /3 




30 


38 


38 


38 


38 




lb 




0.20 


-0.22 


+0.22 


-0.22 


+0.22 




It 




-0.03 





-0.44 





-0.44 




Id 




0.20 


-0.21 


+0.21 


-0.21 


+0.21 




a' 













-2 


-2 


e u 
£d 
\b'\ 
arg(fe') 

a 

b 

\c\ 
arg(c) 

J 33 
V"d 
■^33 
ye 
33 


0.0495(17) 
0.131(7) 
1.04(12) 

90(12)° 
2.17(24) 
1.69(13) 
0.80(16) 
-41(18)° 

0.48(2) 

U. LO 1 ± 1 

0.23(1) 


0.0483(16) 
0.128(7) 
1.07(12) 

91(12)° 
2.27(26) 
1.73(13) 
0.86(17) 
-42(19)° 

0.51(2) 

n 

0.34(2) 


0.0483(18) 
0.102(9) 
1.07(11) 

93(12)° 
2.30(42) 
2.21(18) 
1.09(33) 
-41(14)° 

0.51(2) 

n i&f'Vi 

0.34(2) 


0.0485(17) 
0.127(7) 
1.05(12) 

95(12)° 
2.03(24) 
1.74(10) 
0.81(17) 
-53(10)° 

0.51(2) 

0.34(2) 


0.0485(18) 
0.101(9) 
1.06(10) 

95(12)° 
1.89(35) 
2.26(20) 
1.10(35) 
-41(12)° 

0.51(2) 

n t K&xx\ 

0.34(2) 


(m b /m T )(M x ) 
(3m s /m„){M x ) 
(m d /3m e )(M x ) 


0.67(4) 
0.60(3) 
0.71(7) 


1.00(4) 
0.9(1) 
1.04(8) 


1.00(4) 
0.6(1) 
0.68(6) 


1.00(4) 
0.9(1) 
1.04(8) 


1.00(4) 
0.6(1) 
0.68(6) 




dct Y d (M x ) 
dctV(Mx) 




0.3(1) 


0.92(14) 


0.4(1) 


0.92(14) 


0.4(1) 



Table 4: A x 2 fit of eqs (|15H6[) including the SUSY threshold effects parameterized by the specified 7$. 



phases 4>\ and 02 given by 

0! = (0' 6 - - (<f> b - 0/) (17) 
02 - (0c - 0a) - (<t>b ~ <(>f). (18) 

where <p x is the phase of parameter x. For this reason it is sufficient to consider only b' and c as complex 
with all other parameters real. 

As mentioned above the data favours a texture zero in the (1,1) position. With a symmetric form for the 
mass matrix for the first two families, this leads to the phenomenologically successful Gatto Sartori Tonin 
[54] relation 

V US (M X ) « l&c-l&V^e 

This relation gives an excellent fit to V us with 0i m ± 90°, and to preserve it we take d, d! to be zero in our 
fits. As discussed above, in SU(5) this texture zero leads to the GUT relation Det{M d ) / Det{M l ) = 1 which, 
with threshold corrections, is in good agreement with experiment. In the case that c is small it was shown 
in RRRV that 0i is to a good approximation the CP violating phase 8 in the Wolfcnstcin parameterization. 
A non-zero c is necessary to avoid the relation V u b/V c b — \frn u /m c and with the improvement in the data, 
it is now necessary to have c larger than was found in RRRV 0. As a result the contribution to CP violation 
coming from 02 is at least 30%. The sign ambiguity in 0i gives rise to an ambiguity in c with the positive 
sign corresponding to the larger value of c seen in Tables |3j and HI 

Table [3J shows results from a x 2 fit of eqs (|15|16j) to to the data in Table [2] in the absence of threshold 
corrections. The error, indicated by the term in brackets, represent the widest axis of the la error ellipse 
in parameter space. The fits labeled 'A' have phases such that we have the smaller magnitude solution of 
\c\, and fits labeled 'B' have phases such that we have the larger magnitude solution of \c\. As discussed 
above, it is not possible unambiguously to determine the relative contributions of the off-diagonal elements 
of the up and down Yukawa matrices to the mixing angles. In the fit A2 and B2 we illustrate the uncertainty 

5 As shown in ref. |35l . it is possible, in a basis with large off-diagonal entries, to have an Hermitian pattern with the (1,1) 
and (1,3) zero provided one carefully orchestrates cancelations among Y u and Y d parameters. We find this approach requires 
a strange-quark mass near its upper limit. 



m d , m u, is. 

- W oe 101 

m* V 77v 



(19) 
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associated with this ambiguity, allowing for 0(1) coefficients a'. In all the examples in Tabled the mass 
ratios, and Wolfenstein parameters are essentially the same as in Table [3J 

The effects of the large tan (3 threshold corrections are shown in Table 2J The threshold corrections 
depend on the details of the SUSY spectrum, and we have displayed the effects corresponding to a variety 
of choices for this spectrum. Column A corresponds to a "standard" SUGRA fit - the benchmark Snowmass 
Points and Slopes (SPS) spectra lb of ref([36]). Because the spectra SPS lb has large stop and sbottom 
squark mixing angles, the approximations given in eqns([7l [T0| break down, and the value for the correction 
7i in Column A need to be calculated with the more complete expressions in BRP [IT] . In the column A 
fit and the next two fits in columns B and C, we set a' and c' to zero. Column B corresponds to the fit 
given in the penultimate column of Table [2] which agrees very well with the simple GUT predictions. It is 
characterized by the "anomaly-like" spectrum with M 3 negative. Column C examines the M 3 positive case 
while maintaining the GUT prediction for the third generation m& = m T . It corresponds to the "Just-so" 
Split-Higgs solution. In the fits A, B and C the value of the parameter a is significantly larger than that 
found in RRRV. This causes problems for models based on non-Abelian family symmetries, and it is of 
interest to try to reduce a by allowing a', b' and c' to vary while remaining 0(1) parameters. Doing this for 
the fits B and C leads to the fits B2 and C2 given in Table [4] where it may be seen that the extent to which 
a can be reduced is quite limited. Adjusting to this is a challenge for the broken family-symmetry models. 

Although we have included the finite corrections to match the MSSM theory onto the standard model 
at an effective SUSY scale Ms = 500 GeV, we have not included finite corrections from matching onto a 
specific GUT model. Precise threshold corrections cannot be rigorously calculated without a specific GUT 
model. Here we only estimate the order of magnitude of corrections to the mass relations in Table [2] from 
matching the MSSM values onto a GUT model at the GUT scale. The tan/3 enhanced corrections in eq([7] fTUl) 
arise from soft SUSY breaking interactions and are suppressed by factors of Msusy /Mqut m the high-scale 
matching. Allowing for 0(1) splitting of the mass ratios of the heavy states, one obtains corrections to y b /y T 

2 2 

(likewise for the lighter generations) of 0( u w \z ) from the X and Y gauge bosons and 0( A ) from colored 
Higgs states. Because we have a different Higgs representations for different generations, these threshold 
correction will be different for correcting the Zmsjm^ relation than the mb/m T relation. These factors can 
be enhanced in the case there are multiple Higgs representation. For an SU (5) SUSY GUT these corrections 
are of the order of 2 %. Plank scale suppressed operators can also induce corrections to both the unification 
scale [37] and may have significant effects on the masses of the lighter generations [35J . In the case that the 
Yukawa texture is given by a broken family symmetry in terms of an expansion parameter e, one expects 
model dependent corrections of order e which may be significant. 

In summary, in the light of the significant improvement in the measurement of fermion mass parameters, 
we have analyzed the possibility that the fermion mass structure results from an underlying supersymmetric 
GUT at a very high-scale mirroring the unification found for the gauge couplings. Use of the RG equations 
to continue the mass parameters to the GUT scale shows that, although qualitatively in agreement with 
the GUT predictions coming from simple Higgs structures, there is a small quantitative discrepancy. We 
have shown that these discrepancies may be eliminated by finite radiative threshold corrections involving 
the supersymmetric partners of the Standard-Model states. The required magnitude of these corrections is 
what is expected at large tan /3, and the form needed corresponds to a supersymmetric spectrum in which 
the gluino mass is negative with the opposite sign to the Wino mass. We have also performed a fit to the 
recent data to extract the underlying Yukawa coupling matrices for the quarks and leptons. This is done 
in the basis in which the mass matrices are hierarchical in structure with the off-diagonal elements small 
relative to the appropriate combinations of on-diagonal matrix elements, the basis most likely to be relevant 
if the fermion mass structure is due to a spontaneously broken family symmetry. We have explored the 
effect of SUSY threshold corrections for a variety of SUSY spectra. The resulting structure has significant 
differences from previous fits, and we hope will provide the "data" for developing models of fermion masses 
such as those based on a broken family symmetry. 

M.S. acknowledges support from the United States Air Force Institute of Technology. The views expressed 
in this letter are those of the authors and do not reflect the official policy or position of the United States 
Air Force, Department of Defense, or the US Government. 
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